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We investigate the time evolution of entanglement under various models of decoherence: A general 
heuristic model based on local relaxation and dephasing times, and two microscopic models describ- 
ing decoherence of electron spin qubits in quantum dots due to the hyperfine interaction with the 
nuclei. For each of the decoherence models, we investigate and compare how long the entanglement 
can be detected. We also introduce filtered witness operators, which extend the available detection 
time, and investigate this detection time for various multipartite entangled states. By comparing 
the time required for detection with the time required for generation and manipulation of entangle- 
ment, we estimate for a range of different entangled states how many qubits can be entangled in a 
one-dimensional array of electron spin qubits. 

PACS numbers: 03.65.Ud, 03.67.Mn, 73.21.La 



I. INTRODUCTION 

Entanglement refers to non-classical correlations be- 
tween twol&2& or rnore^ quantum particles, and the 
creation of multiparticle entangled states constitutes a 
key step towards quantum computation 7 -. In this work, 
we investigate the evolution of entangled states under dif- 
ferent models of decoherence: a heuristic model with a 
broad range of applications, and two microscopic models 
specific for electron spin qubits in quantum dots^. We 
show how entanglement can be detected, and how fast 
this needs to be done before the states become disentan- 
gled by decoherence. We also estimate for four different 
classes of multipartite entanglement which class survives 
the longest, and how many entangled qubits can be gen- 
erated and detected with actual experimental means and 
currently known decoherence times. 

We consider decoherence of a local nature, i.e. the 
qubits decohere in an uncorrelated way, as is the case 
in solid-state nanosystems such as electron spin qubits 
in quantum dots^, various superconducting qubits^ and 
other solid-state implementations^: In these systems, 
the decoherence can be characterized phenomenologi- 
cally^i by a relaxation time T\ and a dephasing time T2. 
The microscopic origin of the decoherence is still a mat- 
ter of intensive research: In this paper we discuss some 
microscopic models for electron spin qubits in quantum 
dots and compare them with this heuristic model. 

Our proposed means to analyze entanglement are 
so-called witness operator s 12 ' 13 ! 14 ' 15 ' 16 : locally decom- 
posable observables with a positive expectation value 
for all separable states, and a negative expectation 
value for at least one entangled state. The advantage 
of entanglement witnesses over other methods such as 
e.g. full state tomography is that they require less 
measurements, and thus less experimental effort to 
detect and prove the existence of entanglement for a 
given (mixed) state. Witnesses have intensively been 



used in experiments with photon s 1 and trapped 
ions 1 ^, but so far only few theoretical proposals exist for 
using witness operators in solid-state nanosystems^. 

This paper is organized as follows: in Sec. [II] we briefly 
summarize the mechanisms influencing the time scales T\ 
and T2 for electron spin qubits. We describe two theo- 
retical models of dephasing for these qubits and compare 
them to the heuristic master equation model. We also 
show how we calculate the decoherence of a multipartite 
state of separated qubits using the Lindblad formalism. 
Next (Sec. lIII|) we demonstrate our main ideas and meth- 
ods using the simplest case of two qubits, compare the 
different models of decoherence and introduce our filtered 
witness operator. In Sec. HVI wc consider both GHZ- and 
W-states for three qubits. In Sec. [V] we do the same for 
four qubits and consider in addition the cluster and Dickc 
entanglement classes. We also discuss the dependence of 
the decay of entanglement on the initial state. Finally, 
in Sec. I VII we discuss the case of N qubits: We show 
that the entanglement of a specific GHZ-state can theo- 
retically be detected for any finite time, and discuss the 
feasibility of generating and detecting many-qubit entan- 
gled states of electron spin qubits based on decoherence 
and operation time scales that have been measured in 
recent experiments on single and double quantum dots. 



II. DECOHERENCE MODEL 

Decoherence is caused by uncontrolled interactions be- 
tween the qubit and the environmental. This effect is 
usually characterized by two time scales, the phase ran- 
domization time T2 ( "dephasing time" ) and the time Jf 
in which the excited state |1) relaxes to the ground state 
|0) by energy exchange with the environment ( "relaxation 
time" For electron spin qubits (as for most solid-state 
qubits) the dephasing time is much shorter than the re- 
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laxation time, T2 <C T\, and is therefore the dominant 
time scale for the loss of quantum correlations. In this 
section, we consider both a simple exponential model of 
decoherence based on these two time scales, as well as 
use two microscopic descriptions of dephasing for electron 
spin qubits in quantum dots to derive more sophisticated 
time evolutions of decoherence. 

We start by briefly discussing decoherence mechanisms 
for electron spin qubits. The original idea by Loss and 
DiVincenzo^ proposed to confine single electrons in a 
quantum dot (an island of charge in a two-dimensional 
electron gas) and apply a magnetic field to split the de- 
generacy of the spin-up and spin-down states, thus cre- 
ating a two-level system serving as carrier for quantum 
information: an electron spin qubit. Two electron spin 
qubits interact via a Heisenberg coupling, and this inter- 
action can be controlled by tuning the potential barrier 
between two neighboring dots^ 3 -. Single qubit operations 
rely on electron spin resonance and can be performed by 
applying local electric or magnetic field o 24 ^ 25 . 
Decoherence - interaction with the environment - is 
mainly mediated by two processes, spin-orbit interaction 
and hyperfine interaction^. Spin-orbit interaction does 
not have a direct effect on the electron spin, since the 
electrons do not move, but it leads to a mixing of spin 
and orbital degrees of freedom 2 ^. In GaAs quantum dots, 
spin-orbit interaction is estimated to be small - both ex- 
perimentally 2 ^ and theoretically^ - compared to the hy- 
perfine interaction with the nuclei, so that the latter is 
the dominant source of dephasing. 

If the atoms of the semiconductor material have a non- 
zero nuclear magnetic moment (as for example in GaAs; 
in other materials, such as purified SiC, this effect is not 
present), the electron spin S interacts with the nuclear 
spins via the hyperfine interaction^: the Hamiltonian for 
such a system can be written a a 31 ' 32 

H hf = b S z + e nz I z +h-S. (1) 

Here, b = g* (ibB (e nz = gifJ- n Bo) is the electron (nu- 
clear) Zeeman splitting [calculated using the Bohr (nu- 
clear) magneton /j,b (fJ>N, where /ijv -C /is) and the ef- 
fective (/-factor of the electron (nuclei), g* (gi), which in 
GaAs takes the value g* = -0.44]. Next, I z = YJl=o % 
is the sum over the z-component of all nuclear spins P*, 
and h = X)fc=o denotes the quantum field of the nu- 
clei acting on the electron spin, where n are the number 
of nuclei whose wave function overlaps with the electron's 
wave function (n « 10 5 for typical dots), Ak denotes the 
coupling strength between the fc-th nucleus and the elec- 
tron. Since the electron's wave function is zero outside 
the dot, there is no overlap with the nuclei outside the 
quantum dot - thus each electron in the array couples 
to a different bath of nuclei, and the decoherence is thus 
local, as stated above. Since hyperfine interaction is the 
dominant source of noise, we can disregard other types 
of noise which might induce some correlations between 
different qubits, as for example phonons. 



For an intuitive semiclassical description of decoherence 
due to hyperfine interaction the quantum field of the 
nuclear spins can be treated as an additional (classi- 
cal) magnetic field - the Overhauser field - by replac- 
ing g* [isBn = h. The maximum value this field can 
reach in GaAs is about 3 -^ B max = 5 T for fully polar- 
ized nuclei. In low external magnetic fields, the Over- 
hauser field undergoes Gaussian fluctuations around a 
root-mean-square valu o 31 i 32 ' 33 of B max /y/n. The elec- 
tron thus feels a total magnetic field which consists of the 
sum of the controlled external field Bq and the random 
Overhauser field B n . The field's longitudinal component 
B* uc (parallel to B ) changes the precession frequency 
of the electron spin by h z = g* HbB* iuc . The transverse 
part B*^. changes the precession frequency even only 
in quadratic order, w g* ^i B (B x ' z )^ luc / B . This random 
nuclear magnetic field changes in time: two nuclei with 
different coupling strength Ak can exchange their spin, 
thus leading to a change in the Overhauser field B n ; these 
fluctuations appear on a time scale of 10 — 100 fxs (for a 
weak external field but could probably be extended 
up to well more than several seconds to minutes (for the 
longitudinal nuclear field B^ uc in a strong external field 
Bq)21- 

The bulk dephasing time T 2 (at which the fluctuating nu- 
clear magnetic field removes the phase information) can 
be measured by rotating the spin in the xy-plane, let it 
evolve freely, and then rotate it back along the z-axis for 
measurement (so-called spin echo measurements). Each 
data point then has to be averaged over many measure- 
ments, during which the nuclear field evolves. This leads 
to an average dephasing time T 2 * , which has been mea- 
sured to be about 100 ns 3 ^. 

The dephasing time T 2 of a single electron, on the other 
hand, is very hard to measure, because it is not possi- 
ble to measure the initial orientation and strength of the 
nuclear field with sufficient precision. Estimates in var- 
ious regimes predict T 2 ~ I — 100 pts: a good way to 
estimate T 2 is by using a Hahn echo technique, where 
the free evolution of the spin due to the initial magnetic 
field is undone by reversing the spin, but not the dephas- 
ing due to the change of this nuclear field 4 ^. Assuming 
Gaussian fluctuations of the nuclear field on a time scale 
of 10 s and (a conservatively estimated) T 2 * = 10 ns, a 
time T 2 = 10 /xs can be extracted^, which has been con- 
firmed by measurements^ 3 - providing a lower bound on 
T 2 of 1.2 fjs. 

For a microscopic, quantum-mechanical treatment of 
h, we rewrite Hhf in ([I]) in a parallel and transversal 
part 31 . 32 i 39 ' 40 

H hf = {b + h,)S a + ±(h + S- + h-S+). (2) 



V describes a flip-flop interaction between the electron 
and a nucleus, thus the operators are the raising and low- 
ering operators for the spin (S± — S x ±iS y ) and a nucleus 
(h± = h x ± ihy). This perturbation V is small as soon as 
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there is some external magnetic field and the energy mis- 
match between the electron and the nuclear spin states 
suppresses it, as discussed above: expressed in num- 
bers, this requires \A\ < 2\g*p B B - gipNBo + 2pIA\22-, 
or equivalcntly (using pb Pn and low polarisation) 
|^4| <C 2\g* /j,bBo\, which is fulfilled in typical experi- 
ments (with an external field above > 3 T; in Refsi 35 ' 46 , 
for example, fields of up to 8 T are used). A first ap- 
proximation is to completely neglect this term and only 
consider the change of precession frequency due to the 
nuclear field h z . Using the central limit theorem for a 
large number of nuclear spins results in a Gaussian dis- 
tribution for h z . The transverse correlator, defined as 
the self correlation function of the transverse spin com- 
ponent, (S + ) t = Tr |e iW £/ t S , +e~ iW °'V(0)l (here, p(0) is 
the initial density matrix of the combined system of elec- 
tron and nuclei), is given b y 39 i 40 



(5 (_ se)) 



(5+) exp 



t 2 
2r 2 



■(b+(h z ))t 



(3) 



As opposed to exponential decay of phase coherence with 
time scale T 2 , ^ represents superexponential decay with 
a characteristic time r se = 2H/ ' A^Jnj '(1 — p 2 ): for a GaAs 
dot with almost no polarization, p <C 1, one can estimate 
r sc sa 5 ns, which is much faster than the experimentally 
observed T2-time. The second, imaginary part represents 
the coherent rotation induced by the total magnetic field. 
The value for (h z ) depends on the initial state of the nu- 
clei: for a pure state with each nucleus having probability 
(1 +p)/2 for being in the excited state, it can directly 
be calculated as (h z ) = pA/2, where A is the hypcrfinc 



coupling field. 

A more sophisticated approach is to include the per- 
turbation term V in ((2|), and rewrite the von Neumann 
equation in the form of a Nakajima-Zwanzig generalized 
master equation (GME) 3 ^: 



Pp(t) = -iPLPp(t) - i [ £(i - t')p{t')dt'. 
Jo 



(4) 



Here P is the projector on the electron-subspace, L the 
Liouville-operator (LO = [H, O] for any operator O) and 
£(t) the self-energy superoperator. Using regular per- 
turbation theory in the parameter 1/bo (i.e. for a high 
magnetic field bo 3> A), some (unphysical) secular terms 
arise; these terms do not occur by directly expanding £(i) 
in the GME. The latter results in a self correlation func- 
tion for the transversal spin of the form 3 ^ (in the frame 
oscillating with a frequency proportional to the Zeeman 
splitting) 



{s {bm) } 



(S+) t + R+(t)dt, 



(5) 



where (5+)t is the Markovian solution, and R+(t) is the 
remainder term, i.e. the difference between the expo- 
nential and the non-markovian solution in Born approx- 
imation. This can be written as R+(t) = i J* ^(t — 
t'){S^ m ')t'dt' , and solved by iterating to leading order 
in the parameter 5 = A/ (4N[bo + h z }) (corresponding to 
a high external field, since 6 ~ A/bo). The solution de- 
pends strongly on the wave function of the electron: wc 
assume the electron to have a Gaussian wave function in 
two dimensions, resulting in 



R+(t) ~ -<5(S+) exp 



itAN 



2hT hm (b Q + h z 



t 2 



p sin 

7bm 



Tbr> 



Ip 



TbE 



Tbn 



"Hot 



Here we have defined a characteristic time Tbm = 2nh/A phonons - thus the relaxation can be suppressed by ap- 

(Tbm ~ 1m s f° r GaAs quantum dots). In a realistic set- plying a higher external field. The relaxation time Tihas 

ting, this correction term R+(t) is very small, since S is been measured in experiments to range from 170 ms (at 

very small: in GaAs typically 5 « 1(T 6 . Nonetheless, we B = 1.75 T) to 120 ps (B = 14 T)^&. 

will calculate this correction for completeness. In a phenomenological model of decoherence, the time 

Relaxation of an electron spin qubit is caused by the scales T\ and T 2 are incorporated into a master equation 

same two effects as dephasing: spin-orbit and hypcrfinc model for the density matrix with T\ on the diagonal 

interaction. The required energy for the spin-orbit in- (describing the effect of relaxation) and T 2 on the off- 

teraction to flip the spin of the electron is provided by diagonal (describing dephasing) : 
the phonons in the lattice of the semiconductors forming 

the 2DEG, and can be calculated as a function of the dg 

external magnetic field Bq^. The hyperfine contribution dt 
to relaxation manifests itself as flips of the electron spin 

through exchanging its spin state with a nuclear spin. Qualitatively, the off-diagonal phase components de- 

For increasing external field, the energy mismatch be- crease exponentially with a rate 1/T 2 , and the ground 

tween the nuclear spin states and the electron spin state state Qn becomes populated at the expense of the 

grows, and more and more energy has to be absorbed by excited state Q22 , where the normalization condition 



(1/Tl) Q22 
-(I/T2W 



-(i/r a ) 
-(i/ri) 



Q12 
Q22 



(7) 
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(Tr [p(t)] = 1) has to be fulfilled at any time t © is 
a general phciiomcnological model to describe decohcr- 
ence, and can thus be adjusted to describe decoherence 
for a wide range of systems, but it does not include micro- 
scopic information about the quantum processes causing 
the decoherence. 

We now discuss how to extend these decoherence mod- 
els [Eqs. ([3]), © and (J7J)] to more than one qubit. For 
the exponential decay, this is quite straightforward: we 
rewrite (|7J) in the Lindblad formalism^ using the Lind- 
blad operator L: 



r 2 

+— {2a s ga s - <t s <t s q~- gcr s a s ) . 



(8) 



Here, cr± — 1/2 (a x ± ia y ) and a s = cr_o'+ arc products 
of the Pauli matrices. Comparing the density matrices 
resulting from Eqs. ([7]) and ((5J), we can identify Ti = 1/T\ 
and ri + T-2 = 2/T2. The time evolution of a single qubit 
is then found by solving Cg(t) = dg/dt. To extend (JSJ to 
multipartite states, we write the Lindblad operator for 
the fc-th qubit as C k = l®...®l<g>,C<g>l®...®l, 
where £ is the k-th operator of a total of N. The time 
evolution of the total Af-partite state is then given by 
solving as before CNgif) = dg/dt, with Cn = Ylk=i 
By using this definition we implicitly assumed that the 
decoherence of each qubit is governed by the same Ti and 

r 2 . 

For the two other models, given in ([3]) and we 
construct the density matrix of the entangled states in 
a similar manner. Since the decoherence of the various 
qubits is assumed to be independent, we can just multiply 
the corresponding single matrix entries. For statistically 
distributed quantities (as for example (h z )), we have to 
consider the addition rules for distributions with the cor- 
responding variances, and furthermore we have to take 
into account which contributions have to be conjugated 
(e.g. the precession terms due to the magnetic field). 



III. TWO QUBITS 

Let us first explain our methods and definitions for the 
simple case of two qubits. A separable state g s is defined 
as a state which can be written as a convex combination 
of product statesi*2i4 



(9) 



where g A and g B are states in different subsystems A and 
B and the probabilities pi have to fulfill the normalization 
condition '^2 i Pi = 1. If a state is not of this form, it is 
called entangled. 



We use witness operators 



6.12.13.14.15.16 



to investigate 



the entanglement of various states. An observable W is 
called an entanglement witness if it fulfills the following 
two requirements: 



1. For any separable state g s , the expectation value 
of W is larger than zero: 

Tr [Wq s ] = (W) gs > for all g s separable. 

2. There must be at least one entangled state g e for 
which W has a negative expectation value: 

there exists a g e entangled for which (W) ee < 0. 

Therefore, a measured negative expectation value of the 
witness guarantees that the state is entangled. For the 
experimental implementation, entanglement witnesses 
can be decomposed into local measurements (see also be- 
low), and they usually require much fewer measurements 
than procedures such as full state tomography. Thus 
they are experimentally easier to implement. Finally, it 
should be noted that witnesses can be used to quantify 
entanglement, by giving lower bounds on entanglement 
measures^. 

The witnesses we use in this paper are derived from 
the so-called projector-like witness^ 7 .: 



(10) 



with the constant c standing for the maximum overlap 
between the state \tp) and any separable state. Physi- 
cally, this witness encodes the fact that if a state g has 
a fidelity F = Tr [g |-0) larger than c, then g must be 
entangled. 

We first investigate the time evolution of the Bell 
stated I*") = -^=(|01) - 1 10>). The choice of this Bell- 
state, the singlet state, is motivated by the fact that it is 
the ground state of the quantum system consisting of two 
electron spins in a double quantum dot£, thus it is the 
simplest entangled state that can be created in quantum 
dots. 

The density matrix of a singlet state under exponential 
decay can be found from Eqs. |[7J) and ((8|): 



2[l-a(i)] 

a(t) -P{t) 

-/3(i) a(t) 





(11) 



with the factors a(t) = exp[— t/Ti] = exp[— Tit] for re- 
laxation and (3(t) = exp[-2i/T 2 ] = exp[-(r x + T 2 )t] 
for dephasing. With that, the fidelity 7 F = 
Tr [|*-) (*~| Qv- (t)] is given by 



F(t) = -[a(t)+p(t)\, 



(12) 



and the expectation value of the projective witness for 
l^ - ) is then calculated using (fT0|) . 



Tr [Wsp*- (t)] = (Ws)*- (t) = 5 [1 - a(t) - /3(t)] , 

(13) 

where W5 is the witness for the singlet state, W5 = 1/2— 
l^ - ) (^"l- Figured] shows the decay of entanglement for 
this exponential model of decay of the coherence. 
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FIG. 1: Expectation values for the regular [ (|13p , solid red 
line] and the filtered [ (|20[l , dashed blue line] witness, for 
Ti = 10 3 s _1 and T 2 = 10 6 s _1 . The inset shows a zoom 
into the region where the regular witness becomes positive. 
y is chosen such that the witness expectation value becomes 
minimized at any given time, thus it is time-dependent. 



For the other two models of decoherence [Eqs. © 
and ©], we first have to construct the corresponding 
density matrix: for the relaxation, we keep the exponen- 
tial terms e~*/ Tl , but for the dephasing we use the corre- 
lators presented in the previous section. The first model 
is based on the superexponcntial dephasing from Q for 
each qubit, whose density matrix we label with p se (t). We 
have to calculate the entries |01) and |10), therefore the 
conjugation reverses the phase in ([3]), so the off-diagonal 
dephasing terms [p se ](2,3) = [Pse]* 32 ) (the star * stands 
for complex conjugate) are given in terms of the correla- 
tors (S^) t of the z-th dot (i G {1, 2}) by: 



[Pse](2,3) (t) 



exp 



exp 



' 2 -t 2 



-«/4 1) >-(/4 2) »* 



(14) 



where the second line is for identical statistics of the dots 
(thus with identical characteristic times n = r 2 = r so ). 
Including the boundary condition [/9 se ](2,3)(0) = — 1/2 we 
obtain the density matrix: 



Pseit) 



2[1 - o(t)] 








a(t) 







-Pse{t) 

a(t) 





(15) 



with a(t) = exp[— t/Ti] as before and j3 se {t) = 
[p S e](2,3)(t) from (fT4|) . The witness operator W se for de- 
tecting (|15[) is thus the same as in (|13[) with the replace- 
ment [3(t) — > (3 se . The evolution of the corresponding 
witness is shown in figure [2] a) . 

The third model uses the non-markovian Born approx- 
imation for the decay, Eqs. © and ©. The single elcc- 



i?+(i), where the first exponential term stems from the 
Markovian solution, and the remainder term is given 
in ©. In order to set up the density matrix pbm(^) m the 
non-markovian approximation, we replace [3(t) in pip by 

Ajm(i) = (S+ )t{S+ )*i hi the same way as in the 
super-exponential case. 
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FIG. 2: a) Decay of the entanglement witnesses [without 
(solid) and with (dashed) filter] for the superexponential deco- 
herence model using r sc = 4.4ns (corresponding to p = 0.1). 
The inset shows the advantage of the filter operator, which 
works here as well. 

b) Decay of the entanglement witnesses for the non-markovian 
approximation, including the correction term to the Marko- 
vian solution, for characteristic time Tb m = 1/is, but an en- 
hanced smallness parameter S = 0.1 in order to underline the 
effect of the correction term. 



In the next part we introduce a systematic method to 
enlarge the time interval during which entanglement can 
be detected by the witness operators. This method is 
based on applying local (so-called filtering) operators to 
the witness operator and analyzing the measurement re- 
sults in a different way without requiring more measure- 
ments. Analyzing the witness (|13|) . we see immediately 
that it becomes positive (and hence does not detect the 
entanglement anymore) when after some finite time (3(t) 
becomes smaller than 1 — a(t), though it can be shown 
(by virtue of the PPT-criterior*^, for example) that the 
state Q'Si-it) is entangled for any t < oo, i.e. for any 
/3{t) > 0. 
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Therefore, our goal is to construct a witness opera- 
tor which is able to detect the entanglement in the state 
g<s,-{t) at any time. This can be achieved by a filter 
operator JF, 



(16) 



where the Ti are arbitrary invertible matrices acting on 
individual qubits. Since T is local, application of such 
a filter operator on a state p does not change its en- 
tanglement properties, i.e. qT is entangled, iff g is 
entangled. 

Equivalently, one can apply filter operators to witness 
operators W, and the resulting filtered witness operator 
is then given by 



(17) 



As normalization, we choose Tr [W] = Tr [W r } , to make 
the witnesses' mean values comparable. 

Our goal is now to design a filter Ti such that it in- 
creases the negativity of the witness, i.e. it should in- 
crease the weight of the terms a and (3 in (fT^|) . so that 
the filtered witness can be used to detect entanglement 
during longer times. This can be achieved by the follow- 
ing filter: 



1 

y 



(18) 



with y a positive real number. The normalized filtered 
witness for the singlet state then takes the form 



Wf = 



(fi®f 2 )W s (fi®f 2 ) 

10 

y 2 

y 2 

y 4 



1 



1 



If 



(19) 



and the expectation value is given by 



1 



s /*-(t) 



l + y* 



[1 - a(t) - y 2 p(t)] 



(20) 



Clearly, (Ws)^^(t) ^ s negative if y is chosen large enough 
and time dependent and t < oo (thus f3(t) > 0), and the 
negativity of the witness can be optimized by a suitable 
choice of y = y(t) for a given time t. The remaining 
entanglement [which is not detected by Wg, (fl3|) ] in the 
decohering state can then be detected by measuring this 
filtered witness operator. The effectiveness of the filter 
operator crucially depends on the choice of the singlet 
stated - ) as the initial state: it can easily be shown that 
for the other Bell states |$ ± ) = ^(|00) ± |11», the fil- 
tered witness does not lead to any improvement over the 
regular witness operator. The decay of entanglement in 
our model thus strongly depends on the initial state, even 
within the same basis. 



For the experimental implementation, the witness YVg 
can be decomposed into single-qubit measurements^: 



1 



-[(l + J/ 4 )(l<g>l + cr 2 <g>cr z ) 



(21) 



s 4(1 + y *y 

- (1 - y 4 ){a z <g> 1 + 1 <g) <r z ) + 2y 2 (a x ®a x + a y ® a y )]. 

This decomposition requires three measurement settings 
(namely o~i ® o~i with i £ {x,y,z}) instead of the nine 
settings full state tomography would requir o 14 ! 45 . Similar 
decompositions exist for all other witnesses occurring in 
this pape r 17 i 47 . 

In Figs. [T]and [21 the evolution of the expectation values 
for both the regular witnesses (solid line) and the filtered 
witnesses (dashed line) arc plotted. In the experimentally 
relevant limit (I^ 3> Ti), the advantage of the filter oper- 
ator in an experiment does not manifest itself as strongly 
as would be the case for T2 — I\; however, the princi- 
ple advantage that the entanglement can be detected for 
any finite time is demonstrated in the insets by a zoom 
into the region where the unfiltcred witness becomes pos- 
itive. The filtered witness remains negative, albeit with 
a small, exponentially decaying absolute value, for any 
finite time: this proves that the I 1 ? - ) contains at least 
a very small amount of entanglement at any time under 
our decoherence models; however, it will not lead to a 
significant advantage in an experiment, since the noise 
due to imperfect state preparation and measurement fi- 
delities will render it virtually impossible to measure the 
expectation value with such a high precision. 

From this curve, one can conclude that it becomes dif- 
ficult to detect entanglement after more than a few /is un- 
der realistic conditions assuming any of the decay models 
we have considered (for the superexponential decay even 
after a few ns, but for all models the exact time also de- 
pends on the size of the error bars in a given experiment). 
Consequently, any generation scheme for the Bell states 
which requires a generation time longer than this time 
will probably not work in practice. 



Let us conclude this section by a comparison of the 
three models: The evolution of the expectation values of 
the filtered witnesses is plotted in figure [3] (all in units 
of the respective critical time for the comparison). In 
the non-markovian approximation, we have chosen the 
smallness parameter 6 = 0.1 such that the effect of the 
correction can be shown. This comparison reveals ben- 
efits and drawbacks of each decay model: first, the fil- 
tered entanglement operators are all negative for any fi- 
nite time t < 00. The entanglement in the superexponen- 
tial model decays faster than in the other two, and in the 
non-markovian approximation, we see the effect of the 
power-law tail [the second term in ©] as a periodic re- 
bouncing due to the precession around the z-component 
of the nuclear magnetic field, h z . As a result, the ex- 
pectation value of this witness operator is more negative 
than for the purely exponential decay. In the main plot 
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FIG. 3: Comparison of the evolution of the filtered witness 
operators for exponential decay, superexponential decay, and 
the Born approximation, using 5 = 0.1. The unit of time r is 
the critical time in each model (thus r = T%, r sc , rt, m ). 
For short times, both the superexponential and the non- 
markovian approximation are decaying less strongly than the 
exponential model, but after some time, the witness assum- 
ing exponential decay has a greater negativity. The inset 
compares the longer-time behavior of the three models: the 
non-markovian approximation features a periodic recurrence 
of negative values (due to precession around the nuclear mag- 
netic field) , and the slowest long-time decay (disregarding the 
coherent precession). 



showing the short time evolution scaled by the critical 
time for each model, the differences between the mod- 
els are not as pronounced, and they behave roughly the 
same. 

Let us consider a realistic experiment at this point: In 
an experiment, it is likely that errors will arise due to 
imperfect read-out of the electron spin state a 35 ' 46 , which 
will manifest themselves as error bars on the curves for 
the time evolution. This error will make it unlikely to 
detect the entanglement at longer times in a realistic ex- 
periment. Regarding figure [3] with these errors in mind, 
the evolution within the three models is roughly equiva- 
lent. Another source of possible experimental inaccura- 
cies is the preparation of the initial state - in general it 
will not be the exact aimed-for state, but a mixture of 
states. This mixture will influence the use of the witness 
and the filter operator, which as well strongly depends on 
the nature of the mixture; though the precise influence 
is hard to predict, the filter will always improve the wit- 
ness operator to some degree. However, the creation of 
pure singlet states in double quantum dots has already 
been experimentally achieved in a controllable manner 
and with high probability of success^. Therefore we ex- 
pect that our noiseless results can nevertheless be used to 
give qualitative predictions of the decay of entanglement. 
So far, we have considered two entangled qubits and 
found that their entanglement remains detectable for 
about the same time for three decoherence models. When 
considering generalization to many qubits, we note that 
the exponential model features a big advantage compared 
to the other two, since for this model there is a general 



method for calculating the time evolution of the den- 
sity matrix for an arbitrary number of qubits (see also 
Sec. IVI[) , whereas for the two other models we have to 
construct the density matrix for each new state by hand. 
Therefore, in the following sections, we will use the expo- 
nential model for the generalization to multiple qubits. 



IV. THREE QUBITS 

For three or more particles, the situation is more com- 
plicated, since different classes of multiparticle entangle- 
ment exist 4 ^ 4 ^. 

Let us first discuss the notion of partial separability. 
A state can be partially separable, meaning that some of 
the qubit states are separable, but not all. An example 
for three particles is the state 



(22) 



where \4> AB ) is a (possibly entangled) state of two qubits 
(defined on subsystems A and B), and \4> c ) a state of the 
third qubit (defined on subsystem C). The state \ip bs ) is 
separable with respect to a certain bipartite split, so it is 
called biseparable. A mixed state is biseparable, if it can 
be written as mixture of biseparable pure states. 

If a state is not biseparable, it is genuinely multipar- 
tite entangled. There exist different classes of multipar- 
tite entangled states 4 ^ and the number of entanglement 
classes increases with the number of qubits^. An entan- 
glement class can be defined by the following question: 
given a single copy of two pure states \4>) and \(j>), is it 
possible, at least in principle, to transform into \4>) 
(and vice versa) using local transformations only? Even 
if the probability of success is small? For three qubits, 
for example, two entanglement classes exist, the GHZ 
and the W-class. Every genuine multipartite entangled 
three-qubit state can be transformed into one of the two 
states 4 ^ 



\GHZ 3 ) = -^(|010> + |101» ; 



(23) 



\W 3 ) = ^=(1100) + |010) + |001)), (24) 

but, remarkably, these two states cannot (not even 
stochastically) be transformed into each other, and 
are therefore representatives of different entanglement 
classes. 

Let us now investigate the lifetime of these two states 
using our exponential decoherence model, described be- 
low (jSj) . After calculating the time evolution of the two 
states, we obtain for the corresponding fidelities: 



GHZ 



(*) 



-(exp[-2rii] +cxp[-ri<] 



Fwit) 



+2 exp 
1 



(25) 



(exp[-rii] + 2exp[-(r 1 + T 2 )t]) .(26) 
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From these fidelities, the expectation values of the wit- 
nesses can directly be determined as (WG)p G (t) = 1/2 — 
F G Hz{t) for the GHZ state, and {W w ) PW (t) = 2/3 - 
F w (t) for the W state. 




FIG. 4: Expectation values for the witnesses for the tripartite 
GHZ- and W-states, both regular (solid lines) and filtered 
(dashed), as a function of time t. The projected hlter VV^ P 
(see text) is plotted as well, though it is not better than the 
filtered witness W^. Parameters used are Ti = 10 3 s _1 and 
Y2 = 10 6 s _1 . As in figure [T] y is time dependent and chosen 
so as to minimize the witness expectation value. 



Our next step is to apply the hlter operators to the 
witnesses. This yields the following values for the expec- 
tation value of the filtered witness operators: 



Qw(t) 



2(l + 2y 2 + 2y 4 + y 6 ) 
13 

3 (2 + iy 2 + 6y 4 + 2y 6 ) 



2 + (-3 + 2y 2 ) e- rit + (l - 2y 2 ) e" 2rit - 2y 3 e-^ +r ^ 
2 + (-2 + y 2 )e- r ^ -2y 2 e-^ +r ^ 
I 



(27) 
(28) 



These witness operators can be measured by four (for 
the GHZ state) or five (for the W state) measurement 
settings^, compared to the 27 measurement settings re- 
quired for full state tomography. In principle, the witness 
for the W state can be improved by taking the projec- 
tor onto the subspacc with at most two excitations^, 
(I2 = 1 - I 111) (111 I instead of 1). However, in the 
present case this does not give any improvement, since 
the 1 1 11) { 1 1 1 1 state is not populated. The time evolu- 
tion of the witness expectation values (|27|) and Q28|) arc 
plotted in figure H (for Tj = 10 3 s" 1 and T 2 = 10 6 s" 1 ). 



V. FOUR QUBITS 



The more qubits are added, the more distinct classes 
of entangled states arise. For four qubits, we investigate 



the following four classes 
1 



\GHZ 4 ) 
|C 4 ) 



V2 
1 



(|0101) + |1010)). 



(29) 



(|0101) + |0110) + |1001) - |1010)) , (30) 



(|1000) + 10100) + 10010) + 10001)) , (31) 



V6 



(|0011) + |0101) + |0110) 



|1001) + |1010) + |1100)). 



(32) 



All of these states have been realized in various experi- 
ments for different physical systems^, but so far not in 
solid-state nanosystems. Also, some of their decoherence 
properties have been investigated from different theoret- 
ical perspective o 53 ' 54 . The states \GHZ4) and jWk) are 
the four-qubit versions of the states we have investigated 
for three qubits in the previous section. IC4) is a rep- 
resentative of the so-called cluster class^, important in 
the context of one-way quantum computing^. The Dickc 
stated \Da) is an extension of the W-state and consists 
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of all possible permutations of states containing 2 exci- tations. The fidelity of these states evolves as: 



FaHZi{t) — 


\ (exp[- 


-2rit] 


-1- exp[- 


-2(Ti 


fr 2 )t]), 




Fa{t) = 


\ (exp[- 


-2rit] 


-1- exp[- 


-2(r x 


f r 2 )t] + 


ex P [-(2r 1 + r 2 )i]), 


F Wi (t) = 


\ (exp[- 


-Tit] + 


3 exp[- 


-(r X 4 


-r 2 )t]), 




F D4 {t) = 


\ (exp[- 


-2rit] 


-1- exp[- 


-2(Ti 


f r 2 )t] + 


4exp[-(2r 1 + r 2 )t]) 



(33) 
(34) 
(35) 
(36) 



r 



The corresponding projective witnesses can be found, as 
before, using cl — F(i), with c = 1/2 for the cluster and 
GHZ-states, c = 3/4 for the W-state, and c = 2/3 for the 
Dicke stated. 

Again, filter operations can be applied: the resulting 
formulas are lengthy and therefore not given here. The 
improvement over the regular witness again shows (as 
for the general case of N qubits) that the GHZ-state 
contains in theory entanglement for any finite time - but 
so little, that this result is of a theoretical nature and not 
experimentally relevant. For the other classes of states, 
the filter can lead to a slightly higher negativity, but 
not to an extension of the time where the expectation 
value will become positive. So what is left is to compare 
the differences in the evolution of the expectation values 
of these witness operators for the four classes, and to 
see which one is the most stable, i.e. detectable for the 
longest time. This is done in figure 








- - <W£ 





0.5 



FIG. 5: Expectation values for the filtered witness operators 
for the four considered classes of fourpartite entangled states 
(GHZ, cluster, W and Dicke). The filter operator makes the 
entanglement in the GHZ-state detectable for arbitrarily long 
times (at least in principle). The entanglement of the other 
three classes decays faster, but the cluster state is more sta- 
ble (decays slower) than both the W- and the Dicke states. 
Parameters used are the same as in figure [4] 



At this point the same question can be asked as for 
the two-qubit state that we investigated in Sec. IIIII how 
does the available detection time depend on the exact 
state chosen as representative of a class? Or, equiva- 
lently, the fidelity of which state decoheres most slowly? 
In fact, writing the states above in a different basis leads 
to different decay rates. 

This is illustrated in figure^ where the evolution of the 
witness expectation value of four different cluster states 
is plotted: \C&) from (f30|) . IC^ 16 ' 1 ) is the original clus- 
ter state from Ref^ containing 16 terms, and the two 
additional representations 



IQ 



(|0000) + |0011) + |1110) + 1 1101)) , (37) 



\C. 



(8b 



1 



tfjl (I 0000 ) 



|0011) + |0100) + |0111) 
-| 1000) - |1010) + 1 1101) + |1110)) , (38) 



with 4 and 8 terms, respectively (the first one is a repre- 
sentation with the minimal number of terms, which will 
be used again in the next section, the second one a ro- 
tated version of the original cluster state). As can be seen 
in figure the detection time decreases as the number 
of terms increases, though the effect is not very large for 
four and more terms. The representations with the min- 
imal number of terms thus decohere more slowly. This 
is not surprising, since one can prove for a similar deco- 
herence model that states with the minimal number of 
terms are most robusl^. For representations with the 
same number of terms, the number of excitations in each 
term can influence the detectability: which one of the 
two is easier to detect then depends on the ratio of Ti 
and r 2 . 



VI. N QUBITS 

Let us now consider the general situation of TV qubits. 
We concentrate on three types of entangled states for 
which there exist proposals how to generate them using 
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,56 



— c 4 


■ cf 


■cf 


- - cf 



\C N ) (Cat I 



(44) 



fe=i 



with Sk a product of Pauli matrices. We incorporate the 
effects of dephasing (disregarding the relaxation of the 
qubits, i.e. setting Ti = 0, which leads to an error of 
less than 0.01%o for four qubits) in every term in the 
sum of the expanded (14"4")) . The resulting fidelity of the 
cluster state can then easily be calculated numerically up 
to N = 24 qubits. 



FIG. 6: Expectation values of the witnesses for some four- 
partite cluster states. The upper index indicates the number 
of terms in the representation of each state, the state without 
index is the one given in (|30|l . Parameters used are the same 
as in figure [4] 



available single- and two-qubit operations in quantum 
dots^S: GHZ-, W-, and cluster states. Our goal is to 
calculate the time evolution of the normal (unfiltered) 
witness for arbitrary N and compare this with the time 
necessary to generate and measure the state. The rep- 
resentatives of the first two classes can be written down 
straightforwardly (for even N): 



\GHZ 



N 



1 

7T 



|01... 01) + |10...10)), 



(39) 



\W N ) = -= (|00 . . . 01) + |00 . . . 10) + . . . + 1 10 . . . 00)) . 
v N 

(40) 

Calculating the expectation value of the witness opera- 
tors leads to (for even N): 



(W G ) eG(t) =- <j 1 



-rit 



exp 



N 



(Ww) gw (t) =-{iV-l-exp[-ri*] 

-(iv-i)exph(r 1 + r 2 )t]} 



(Ti + r 3 )t 

(41) 
(42) 



The general form of the cluster state - the one we 
consider here containing the minimal number of terms, 
namely 2 N / 2 - is more complicated; it can be written as^ 



\C 



N 



fc=l 



[|OO) + |ll)(q :c 0l)] 
V2 



(43) 



where this formula should be understood as an iteration, 
with the operator (a x ® 1) acting on the Bell state of 
the next two qubits. For four qubits, this results exactly 
in the representation |C^ 4 ^ ) from (|3"T|) . the evolution of 
which is plotted in figure [5] To calculate the time evo- 
lution of the fidelity, we represent the cluster state (14"3"1) 




FIG. 7: The dots show the time at which the expectation 
value of the (unfiltered) witness for a given class of states 
becomes positive as a function of the qubit number N. Pa- 
rameters used are I\ = 10 3 s _1 and Y2 = 10 6 s _1 , except for 
the cluster state where the relaxation is disregarded (ri = 0). 
Also plotted are (shown by lines) the times required to gen- 
erate the given entangled states, see the text for explanation. 



Figure [7] shows the time at which the expectation value 
of the (unfiltered) projective witness for each of the three 
states [Eqs. (|4ip . (|4"2"|) and (|44p] becomes positive as a 
function of the qubit number N, as well as a rough esti- 
mate of the time necessary to generate and measure these 
states. For electron spin qubits in quantum dots the gen- 
eration times are taken from Ref*^: for both the clus- 
ter and the W-states the time required to produce these 
states is independent of the number of qubits, whereas 
the production time of states of the GHZ-class scales lin- 
early with the number of qubits. The measurement times 
are composed as follows: measurement distinguishes be- 
tween spin-up and spin-down (defined along the z-axis^) 
and measuring the components a x and a y then requires 
a rotation of the spins by 7r/2, which takes about ~ 50 
ns^. The sum of the generation and the measurement 
time is given by the lines in the plot. 

We see in figure [7| that (as in Figs. 5] and [5]) the entan- 
glement of the GHZ-state can be detected for the longest 
times, but it is more time-consuming to generate than the 
other two entangled states. Based on the estimates in fig- 
ure generation and detection of GHZ states should be 
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possible for up to 14 qubits (with the standard projective 
witness and assuming current operation and decoherence 
times for electron spin qubits). The cluster state is the 
state which can be detected for the largest number of 
qubits, although for up to 12 qubits the "time reserve" 
(i.e. the difference between the time needed for genera- 
tion and measurement and the time when the expectation 
value of the witness operator becomes positive) for the 
GHZ state is somewhat larger than for the cluster state. 
The W-state is the least suitable, the largest state would 
contain about ~ 10 qubits. 

Our results for the cluster state show that one-way 
quantum computing^ is not really feasible in quantum 
dots with current dephasing times: we expect that up 
to maximally ~ 12 qubits could be entangled under the 
presented preparation scheme, which is far too few for 
exploiting the advantages of a quantum computer. 

Based on our assumptions, thus, the simplest state to 
generate and prove it's entanglement would be the GHZ- 
state for up to twelve qubits, and the cluster state for 
more than twelve and up to twenty qubits, though the 
remaining entanglement becomes very small. The same 
holds for the filtered witness for the GHZ state for an 
arbitrary number of qubits. 



which locally undergo decoherence. We have considered 
three different phcnomcnological models for the dephas- 
ing of the qubits based on exponential and supcrcxpo- 
nential decay. Using witness operators as detectors of 
entanglement and introducing a specific class of filtered 
witness operators, we estimated the maximum available 
detection time for entanglement of two electrons using 
each of the models and found that the time during which 
entanglement is detectable is independent of the model 
chosen. We then expanded the exponential model to the 
case of multipartite entanglement: For three and four 
qubits, we compared the decay of entanglement for dif- 
ferent classes of entangled states with each other, namely 
the GHZ-, W-, cluster and Dicke classes. We also gave 
limits on the maximum number of entangled qubits that 
can be created and measured based on currently known 
decoherence times for electron spin qubits. The most 
suitable entangled state turns out to be the GHZ-state 
for up to a few qubits. Our results can help to make 
a choice as to which state to prepare in experiments. 
Since local decoherence is characteristic for many types 
of solid-state qubits, our model and the filtered operator 
technique are applicable to a variety of these qubits. 



VII. CONCLUSION 

This work has been supported by The Netherlands 
In conclusion, we have investigated entanglement and Organisation for Scientific Research (NWO), the FWF 
its detectability in a linear array of electron spin qubits (START prize) and the EU (SCALA, OLAQUI, QICS). 



1 A. Einstein, B. Podolski, and N. Rosen, Phys. Rev. 47, 
777 (1935). 

2 E. Schrodinger, Naturwissenschaften 23, 807 (1935). 

3 J. S. Bell, Physics 1, 195 (1964). 

4 R. F. Werner, Phys. Rev. A 40, 4277 (1989). 

5 R. Horodecki, P. Horodecki, M. Horodecki, and K. 
Horodecki, |arXiv:quant-ph/ 0702225 

6 M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. 
A 223, 1 (1996). 

7 M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information, (Cambridge University Press, 
Cambridge, UK, 2000). 

8 R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 
and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 1217 
(2007). 

9 G. Wendin and V. S. Shumeiko, in the Handbook of Theo- 
retical and Computational Nanotechnology, 1, 129; see also 
|arXiv:cond-mat/0508729| 

10 For a recent overview of more solid-state qubit systems 
such as e.g. nitrogen-vacancy centers in diamond, see R. 
Hanson and D. D. Awschalom, Nature 453, 1043 (2008) 
and references therein. 

11 See e.g. H.-A. Engel and D. Loss, Phys. Rev. B 65, 195321 
(2002). 

12 B. M. Terhal, Phys. Lett. A 271, 319 (2000). 

13 M. Lewenstein, B. Kraus, J. I. Cirac, and P. Horodecki, 
Phys. Rev. A 62, 052310 (2000); B. M. Terhal, Phys. Lett. 



A 271, 319 (2000). 

14 O. Giihne, P. Hyllus, D. Brufi, A. Ekert, M. Lewenstein, 
C. Macchiavello, and A. Sanpera, Phys. Rev. A 66, 062305 

(2002) . 

15 O. Giihne, P. Hyllus, D. Bruss, A. Ekert, M. Lewenstein, 
C. Macchiavello, and A. Sanpera, J. Mod. Opt. 50, 1079 

(2003) . 

16 O. Giihne and G. Toth, Physics Reports 474, 1 (2009). 

17 M. Bourennane, M. Eibl, C. Kurtsiefer, S. Gaertner, H. 
Weinfurter, O. Giihne, P. Hyllus, D. Brufi, M. Lewenstein, 
and A. Sanpera, Phys. Rev. Lett. 92, 087902 (2004). 

18 C.-Y. Lu, X.-Q. Zhou, O. Giihne, W.-B. Gao, J. Zhang, 
Z.-S. Yuan, A. Goebel, T. Yang, and J.-W. Pan, Nature 
Physics 3, 91 (2007). 

19 D. Leibfried et ai, Nature 438, 639 (2005); H. Haffner et 
al, Nature 438, 643 (2005). 

20 M. Blaauboer and D. P. DiVincenzo , Phys. Rev. Lett. 95, 
160402 (2005); L. Faoro and F. Taddei, Phys. Rev. B 75, 
165327 (2007). 

21 M. Schlosshauer, Rev. Mod. Phys. 76, 1267 (2005). 

22 D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

23 JR. Petta et ai, Science 309, 2180 (2005). 

24 F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, 
K. C. Nowack, T. Meunier, L. P. Kouwenhoven, and 
L. M. K. Vandersypen, Nature 442, 766 (2006). 

25 K.C. Nowack et ai, Science 318, 1430 (2007). 



12 



26 A. V. Khaetskii and Yu. V. Nazarov, Phys. Rev. B 64, 
125316 (2001). 

27 D. M. Zumbhl, J. B. Miller, C. M. Marcus, K. Campman 
and A. C. Gossard, Phys. Rev. Lett. 89, 276803 (2002), 
J. B. Miller, D. M. Zumbhl, C. M. Marcus, Y. B. Lyanda- 
Geller, G. Goldhaber-Gordon, K. Campman and A. C. 
Gossard, Phys. Rev. Lett. 90, 076807 (2003). 

28 V. N. Golovach, A. Khaetskii, and D. Loss, Phys. Rev. 
Lett. 93, 016601 (2004). 

29 A. Abragam, The Principles of Nuclear Magnetism, (Ox- 
ford University Press, 1961). 

30 D. Paget, G. Lampel, B. Sapoval, and V. I. Safarov, Phys. 
Rev. B 15, 5780 (1977). 

31 A. V. Khaetskii, D. Loss, and L. Glazman, Phys. Rev. Lett. 
88, 186802 (2002). 

32 I. A. Merkulov, A. L. Efros, and M. Rosen, Phys. Rev. B 
65, 205309 (2002); 

33 P.-F. Braun et al, Phys. Rev. Lett. 94, 116601 (2005); 
A. C. Johnson, J. R. Petta, J. M. Taylor, A. Yacoby, M. D. 
Lukin, C. M. Marcus, M. P. Hanson, and A. C. Gossard, 
Nature 435, 925 (2005); F. H. L. Koppens, J. A. Folk, J. M. 
Elzerman, R. Hanson, L. H. Willems van Beveren, I. T. 
Vink, H. P. Tranitz, W. Wegscheider, L. P. Kouwenhoven, 
L. M. K. Vandersypen, Science 309, 1346 (2005); J.M. 
Taylor et al, Phys. Rev. B 76, 035315 (2007). 

34 R. Hanson, B. Witkamp, L. M. K. Vandersypen, L. H. 
Willems van Beveren, J. M. Elzerman, and L. P. Kouwen- 
hoven, Phys. Rev. Lett. 91, 196802 (2003). 

35 J. M. Elzerman, R. Hanson, L. H. Willems van Beveren, B. 
Witkamp, L. M. K. Vandersypen, and L. P. Kouwenhoven, 
Nature 430, 431 (2004). 

36 R. G. Shulman, B. J. Wyluda, and H. J. Hrostowski, Phys. 
Rev. 109, 808 (1958). 

37 This is not confirmed experimentally, but the very slow de- 
cay of nuclear spin polarization of up to several minutes is 
an indication for slow change of B^ uc ; see e.g. A. K. Hiittel, 
J. Weber, A. W. Holleitner, D. Weinmann, K. Eberl, and 
R. H. Blick, Phys. Rev. B 69, 073302 (2004). 

38 J. M. Kikkawa and D. D. Awschalom, Phys. Rev. Lett. 80, 
4313 (1998). 

39 W.A. Coish and D. Loss, Phys. Rev. B 70, 195340 (2004). 

40 L. Chirolli and G. Burkard, Advances in Physics 57, 225 
(2008), see also larXiv:0809.4716l . 

41 B. Herzog and E. L. Hahn, Phys. Rev. 103, 148 (1956). 

42 G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 

43 F. G. S. L. Brandao, Phys. Rev. A 72, 022310 (2005); D. 
Cavalcanti and M. O. Terra Cunha, Appl. Phys. Lett. 89, 
084102 (2006); O. Giihne, M. Reimpell, and R. F. Werner, 
Phys. Rev. Lett. 98, 110502 (2007); J. Eisert, F. Brandao, 
and K. Audenaert, New J. Phys. 9, 46 (2007). 

44 A. Peres, Phys. Rev. Lett. 77, 1413 (1996); M. Horodecki, 
P. Horodecki, and R. Horodecki, Physics Letters A 223, 1 
(1996). 

45 G. Toth and O. Giihne, Phys. Rev. Lett. 94, 060501 



(2005). 

46 R. Hanson, L. H. Willems van Beveren, I. T. Vink, J. M. 
Elzerman, W. J. M. Naber, F. H. L. Koppens, L. P. 
Kouwenhoven, and L. M. K. Vandersypen, Phys. Rev. Lett. 
94, 196802 (2005). 

47 O. Giihne and P. Hyllus, Int. J. Theor. Phys. 42, 1001 

(2003) . 

48 W. Diir, G. Vidal, and J. I. Cirac, Phys. Rev. A 62, 062314 
(2000). 

49 A. Acm, D. Brufi, M. Lewenstein, and A. Sanpera, Phys. 
Rev. Lett. 87, 040401 (2001). 

50 Daniel M. Greenberger, Michael A. Home, Anton 
Zeilinger: Bell's theorem, Quantum Theory, and Concep- 
tions of the Universe, pp. 73-76, (Kluwer Academics, Dor- 
drecht, The Netherlands, 1989), see also larXiv:0712.0921l 
D. M. Greenberger, M. A. Horne, A. Shimony, and A. 
Zeilinger, Amer. J. Phys. 58, 1131-43 (1990). 

51 F. Verstraete, J. Dehaene, B. De Moor, and H. Verschelde, 
Phys. Rev. A 65, 052112 (2002); L. Lamata, J. Leon, D. 
Salgado, and E. Solano, Phys. Rev. A 75, 022318 (2007). 

52 See e.g. Refs.- 1 ^^ andi 9 - for realizations of GHZ, W, and 
cluster states of up to 8 qubits using photons and ions, 
respectively, and N. Kiesel, C. Schmid, G. Toth, E. Solano, 
and H. Weinfurter, Phys. Rev. Lett. 98, 063604 (2007) for 
a four-qubit Dicke photon state. 

53 O. Giihne, F. Bodoky, and M. Blaauboer, Phys. Rev. A 
78, 060301(R) (2008)'. 

54 C. Simon and J. Kempe, Phys. Rev. A 65, 052327 (2002); 
W. Diir and H.J. Briegel, Phys. Rev. Lett. 92, 180403 

(2004) ; L. Aolita, R. Chaves, D. Cavalcanti, A. Acin, and 
L. Davidovich, Phys. Rev. Lett. 100, 080501 (2008); A. 
Borras, A. P. Majtey, A. R. Plastino, M. Casas, and A. 
Plastino, Phys. Rev. A 79, 022108 (2009). 

55 R. Raussendorf and H.-J. Briegel, Phys. Rev. Lett. 86, 
5188 (2001), R. Raussendorf, D. E. Browne and H.-J. 
Briegel, Phys. Rev. A 68, 022312 (2003). 

56 M. Hein, W. Diir, J. Eisert, R. Raussendorf, M. Van den 
Nest, and H.-J. Briegel: in Proceedings of the International 
School of Physics Enrico Fermi on Quantum Computers, 
Algorithms and Chaos, (Varenna, Italy, 2005); see also 
|arXiv:quant-ph/0602096| 

57 R. H. Dicke, Phys. Rev. 93, 99 (1954). 

58 G. Toth, J. Opt. Soc. Am. B 24, 275 (2007). 

59 See (2) in H.-J. Briegel and R. Raussendorf, Phys. Rev. 
Lett. 86, 910 (2001). 

60 F. Bodoky and M. Blaauboer, Phys. Rev. A 76, 052309 
(2007). 

61 In general, a cluster state is defined as an eigenstate of 
certain local observable a 55 ' 56 ; here we choose them as 
Si = a z a z l...; 5 2 = cr x a x cr x l . . .; S3 = \a z a z a z \ . . .; 
. . .; Sn — 1 ■ ■ ■ la x a x . This is locally equivalent to the def- 
inition as i n 55 ' 56 ; however, the state given in (|43[) contains 
the minimal number of terms; see Ref.— . 



